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Abstract 

Selected  results  from  an  extensive  oscillating- 
airfoil  experiment  are  analyzed  and  reviewed.  Four 
distinct  regimes  of  viscous-inviscid  interaction 
are  identified,  corresponding  to  varying  degrees  of 
unsteady  flow  separation.  The  dominant  fluid 
dynamic  phenomena  are  described  for  each  regime. 

Ten  specific  test  cases,  including  the  appropriate 
flow  conditions  and  experimental  results,  are  pro¬ 
posed  for  evaluating  unsteady  viscous  theories  and 
computational  methods. 


ntroduction 


The  free-stream  Mach  number  was  varied  between 
0.07  and  0.30,  with  important  consequencies; 1  how¬ 
ever,  only  the  data  at  H,  '  0.30  are  presented 
here.  For  this  Mach  number,  the  Reynolds  number  was 
approximately  4  x  106  based  on  chord,  and  the 
reduced  frequency  parameter,  k  =  wc/2U„,  varied  from 
0  to  0.20. 

Details  of  the  instrumentation,  experimental 
procedures,  measurement  uncertainties,  and  test  con¬ 
ditions,  as  well  as  the  coordinates  of  each  airfoil, 
are  documented  in  Ref.  2.  That  reference  also 
describes  the  computer  data  tapes,  which  will  be 
made  available  upon  request. 


Important  viscous  phenomena  combine  with 
unsteady  effects  in  a  variety  of  current  aeronauti¬ 
cal,  hydrodynamic,  and  wind-energy  problems.  Con¬ 
sequently,  several  recent  investigations  have 
attempted  to  predict,  compute,  or  empirically  cor¬ 
relate  the  unsteady  airloads  on  oscillating  airfoils 
that  experience  varying  degrees  of  flow  separation, 
and  additional  work  is  in  progress  at  a  number  of 
research  centers.  The  results  are  mixed,  and  con¬ 
flicting  claims  abound;  this  is  partly  because  the 
difficulty  of  the  problem  increases  rapidly  with 
increasing  amounts  of  separation,  and  partly  because 
definitive  data  over  wide  ranges  of  separated  flow 
conditions  do  not  exist. 

The  authors  and  their  colleagues  have  recently 
concluded  an  extensive  experimental  investigation  of 
boundary-layer  transition,  separation,  and  unsteady 
stall  on  oscillating  airfoils  in  two-dimensional 
subsonic  flow.  The  experiment  included  more  than 
50  combinations  of  and  parameters  of  the 

unsteady  motion  for  each  of  eight  different  airfoil 
sections.  Highlights  of  the  principal  results, 
insofar  as  the  extensive  lift,  drag,  and  pitching 
moment  data  are  concerned,  were  described  recently 
in  Ref.  1.  The  purpose  of  the  present  paper  is  to 
extract  from  a  limited  number  of  these  test  points, 
the  specific  details  and  experimental  information 
that  might  serve  to  guide  the  development  of  new 
unsteady  viscous  theories  and  computational  methods. 
A  hierarchy  of  model  problems  can  be  developed, 
corresponding  to  Increasing  amounts  of  interaction 
between  the  viscous  and  inviscid  parts  of  the  flow 
field.  Within  each  level,  specific  sets  of  data  are 
proposed  that  can  be  used  as  test  cases  for  evaluat¬ 
ing  existing  and  future  prediction  methods. 

The  data  for  the  present  paper  were  obtained 
from  the  three  airfoils  shown  in  Fig.  1.  Pressure 
and  hot-wire  instrumented  models,  0.62-m  chord  and 
2.1-m  span,  were  oscillated  in  pitch, 
i  =  i(J  +  ,,  sin  mt,  about  an  axis  at  x/c  =  0.25 
in  a  2.1-  by  3.0-m  atmospheric-pressure,  solid-wall 
wind  tunneL.  The  effects  of  the  tunnel  walls  are 
thought  to  be  small,  but  not  negligible. 
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Summary  of  the  Viscous-inviscid 
Interaction  Regimes 


The  overall  features  of  the  flow  field  around 
an  airfoil  In  subsonic  flow  are  primarily  charac¬ 
terized  by  the  degree  or  extent  of  flow  separation. 
For  a  given  airfoil,  the  primary  parameter  that 
determines  the  degree  of  separation  is  the  maximum 
angle  of  attack  (amax  =  a0  +  .  i  for  sinusoidal 
oscillations).  An  important  aspect  of  the  flows 
discussed  in  this  paper  is  the  large  amplitudes,  ij 
that  produce  the  large  maximum  angles.  This  con¬ 
trasts  with  the  hierarchy  of  viscous  effects  on 
oscillating  airfoils  at  transonic  speeds  and  low 
angles  of  attack, 3,‘*  where  the  scale  of  the  inter¬ 
action  is  governed  primarily  by  the  strength  and 
motion  of  the  shock  wave.  So  far,  prediction 
methods  for  this  class  of  problems  have  not  been 
successful  for  the  low-speed,  high-angle  problems, 
and  vice  versa. 


The  importance  of  iq^x  is  illustrated  in 
Fig.  2,  which  portrays  four  important  regimes  of 
viscous-inviscid  interaction  for  oscillating  air¬ 
foils.  For  the  left-hand  part  of  the  figure,  that 
is,  =  13°,  there  was  almost  no  separation 

throughout  the  cycle,  although  unsteady  boundary- 
layer  displacement  thickness  effects  were  not  com¬ 
pletely  negligible.  When  cqxax  was  increased  to 
14°,  the  limited  separation  that  occurred  during  a 
small  fraction  of  the  cycle  distorted  the  hysteresis 
loops  of  the  unsteady  pressures  and  airloads,  from 
a  practical  standpoint,  the  effect  on  is 

particularly  important.  This  stall-onset  condition 
represents  the  limiting  case  of  the  maximum  lilt 
that  can  be  obtained  with  no  significant  penalty  in 
pitching  moment  or  drag. 


A  slight  additional  Increase  in  to  1  3" 

produced  a  major  increase  in  the  extent,  severity, 
and  duration  of  the  separation  phenomenon,  for  the 
conditions  shown  in  Fig.  2.  This  type  of  viscous- 
inviscid  interaction  produced  what  is  called  light 
dynamic  stall.1  Further  increases  in  imilx  led  to 
the  deep  dynamic  stall  regime,  with  a  large  viscous 
zone  over  the  entire  upper  surface  of  the  airfoil 
during  half  or  more  of  the  cycle. 

Within  each  of  the  stal  I  regimes  shown  in 
Fig.  2,  the  details  of  the  flow  field  depend  on  Mach 
number,  Reynolds  number,  airfoil  shape,  and  the 


I 


airfoil  motion,  a(t).  To  a  first  approximation,  M„ 
and  leading-edge  geometry  are  the  principal  deter¬ 
minants  of  the  type  of  boundary-layer  separation 
(for  example,  leading-edge,  trailing-edge,  or  mixed 
separation  characteristics),  although  changes  from 
one  type  to  another  can  occur  as  the  frequency  of 
the  oscillation  increases.1  For  a  given  class  of 
separation  or  stall  behavior,  the  amplitude  and 

reduced  frequency  k  «  uic/2U„  largely  determine  the 
stall  and  reattachment  angles  and  the  size  and  shape 
of  the  hysteresis  loops;  or  alternatively,  the  depar¬ 
ture  from  quasi-static  behavior. 

More  detailed  descriptions  of  the  various  stall 
mechanisms  can  be  found  in  Ref.  1,  Refs.  5-7,  and 
elsewhere.  However,  it  is  apparent  from  the  brief 
remarks  above  that  passing  from  no  stall  to  deep 
dynamic  stall  encompasses  a  wide  range  of  viscous 
flow  phenomena,  and  that  realistic  but  efficient 
prediction  methods  should  probably  be  tailored  to 
the  specific  characteristics  of  the  flow  regime  of 
interest.  The  following  sections  provide  some  spe¬ 
cific  experimental  Information  and  test  cases  that 
can  be  used  in  developing  and  validating  new  anal¬ 
yses  and  computations  for  each  regime  depicted  in 
Kig.  2. 


Ill,  No  Stall;  Weak  Interactions 
General  Features 

If  viscous  effects  are  confined  to  thin  bound - 
a-y  layers,  the  primary  effects  of  oscillation  can 
be  derived  from  unsteady  thin-airfoil  theory.  The 
first-order  results  for  pressure  and  lift  take  the 
following  form: 

cp  y— fi(x’ c)  « 


C[  =  --.zzzrzrz  f,(l)  (2) 

.1  -  M‘ 

where  the  minus  and  plus  signs  refer  to  the  upper 
and  lower  surfaces,  respectively.  All  of  the 
unsteady  effects,  such  as  amplitude  and  phase  changes 
with  respect  to  the  motion  u(t),  are  contained  in 
the  functions  ft  and  f , .  The  remaining  terms  in 
Kqs.  (I)  and  (2)  are  easily  recognized  as  the  solu- 
lion  for  a  flat  plate  in  steady  flow,  and 
fj  *  f>  ►  1  in  the  quasi-steady  limit.  The  details 
of  these  unsteady  solutions  are  readlty  available  in 
the  vast  literature  on  thin-airfoil  theory  (e.g.. 

Kef.  8). 

A  number  of  secondary  effects,  such  as  airfoil 
thickness  and  camber,  finite  mean  angle  of  attack, 
large-amplitude  motion,  boundary- layer  displacement 
thickness  effects  (including  small  amounts  of  separa¬ 
tion),  and  wind-tunnel  wall  corrections,  can  be 
obtained,  at  least  qualitatively,  by  superposition 
of  individual  corrections  to  the  unsteady  flat-plate 
solution.  This  is  indicated  schematically  in  Fig.  3 
for  the  surface  velocity,  U.  For  quantitative  pre¬ 
cision,  the  pressure  should  be  calculated  from  U 
and  using  the  nonlinear  Bernoulli  equation. 

Although  conceptually  rather  simple,  accurate  pre¬ 
dictions  of  these  effects  are  not  trivial.  Conse¬ 
quent  tv,  experimental  test  cases  in  this  weak  inter¬ 
act  Ion  regime  represent  the  first  level  in  the 
hierarchy  of  model  problems  to  be  considered.  The 


following  four  examples  are  arranged  in  the  order  of 
increasing  levels  of  complexity. 

Case  1:  NACA  0012,  u  ■  5°  +  5°  sin  u>t 

Figure  4  illustrates  the  basic  features  of  an 
oscillating  airfoil  with  omax  rather  large  but 
well  below  the  static  stall  angle.  The  boundary 
layers  on  both  the  upper  and  lower  surfaces  were 
fully  attached  throughout  the  cycle,  except  for  a 
small  separation  bubble  near  the  upper-surlace 
leading  edge  for  a  i  5°,  which  produced  transition 
from  laminar  to  turbulent  flow. 

Figure  5  shows  the  details  of  the  unsteady 
pressure  distributions  by  harmonic  components.  The 
Cp  values  are  plotted  vs  /x  rather  than  x,  so  as 
to  stretch  out  the  leading-edge  region  where  the 
variations,  and  the  discrepancies  with  linear 
theory,  are  the  greatest.  Linear  theory  predicts 
C  “  0  and  equal  and  opposite  values  for  the 
individual  components  of  Cp0  and  Cpi  on  the  upper 
and  lower  surfaces,  which  is  clearly  not  the  case 
in  reality.  However,  the  general  trends  of  the 
measured  pressures  are  precisely  those  suggested  by 
the  linear  superposition  concepts  indicated  in 
Fig.  3,  provided  the  nonlinear  Bernoulli  equation 
is  used,  as  mentioned  above. 

It  is  interesting  to  note  that  linear  iheorv 
gives  a  much  better  estimate  of  the  r 

between  the  upper  and  lower  surface  pressures,  than 
it  does  of  the  single-surface  values.  As  a  result, 
the  mean  and  fluctuating  values  of  the  lilt  coeffi¬ 
cient  given  by  Eq.  (2)  are  approximately  correct 
for  this  case. 

Finally,  there  is  a  point  to  he  mad*'  about  a 
relatively  minor  aspect  of  the  lift  behavior  shown 
in  Fig.  4.  The  hysteresis  in  Cj  versus  *  is 
slightly  different  from  that  predicted  by  l tin  at 
theory.  This  small  difference,  also  evident  in  tin 
results  on  the  other  profiles,  might  be  as,  i  il-.d  t, 
the  effect  of  airfoil  thickness,  but  it  is  ai-,,-  i  h. 
trend  predicted  for  wall  interference  %-tt.it-..  ‘ 
Computations  or  analyses  of  all  of  tin  .  as.  •.  m  :‘n 
paper  should  Include  the  solid  wind  tuini.  l  wall- 
outer  boundaries,  rat  her  than  t  i  .» i  i  ..-n.iiii'n 


Case  2:  M.K-7301 , 


♦  sin  i 


Pilfering  from  Use  l  onlv  in  i 
Case  2  produced  the  data  plot  lt  d  m 
effects  of  both  thickness  and  «  ambei 
in  the  mean  pressure  distribution,  , 
expected.  However,  the  fluctuating 
to  be  much  less  sensitive  to  the  !u 
and  are  mainly  affected  by  t lie  Mum 
geometry.  The  leading-edge  radius  i 
of  the  0012  airfoil,  and,  consequent 
values  of  Cpi  ne.u  the  leading  edg 
NLR-7301  airfoil.  It  should  be  note 
theory  for  this  case  is  identical  to 
Case  1 . 


Case  3:  NACA  0012, 


8 "  t  S  1  sin 


Figure  7  shows  the  principal  results  |.„  ihi-. 
more  complex  ruse.  Increasing  <*,.,*  h.is  two  main 
effects,  compared  with  Case  l.  first,  the  Ou,  tun¬ 
ing  pressures  Increase  due  to  the  non II neat  o.iilti 
but  Ions  of  mean  angle  effects.  Second,  the  upp.  t 
surface  boundary  layer  thickens,  especially  neat  t he 
trailing  edge.  The  former  el  led  is  most  pronoumed 


« 


for  the  in-phase  pressure  distribution  near  the 
leading  edge,  whereas  the  latter  is  the  main  reason 
for  the  change  in  the  out-of -phase  component  of  Cp. 
This  change  in  quadrature  pressure  is  also  reflected 
in  the  sense  of  the  -  a  hysteresis  loop  (not 

shown) ,  which  was  opposite  to  that  of  the  two  pre¬ 
vious  cases.  Again,  the  linear-theory  curves  of 
Cp/u,  omitted  for  clarity  in  Fig.  7,  would  be  iden¬ 
tical  to  those  shown  in  Fig.  5  for  Case  1. 

Case  4:  NLR-7301,  a  -  10°  +  5°  sin  u>t 

In  Case  A,  a  thin  layer  of  reversed  boundary- 
layer  flow  was  detected  moving  forward  from  the 
trailing  edge  to  approximately  midchord  as  a 
approached  <%££•  The  front  edge  of  this  reversed- 
flow  layer  then  retreated  toward  the  trailing  edge 
as  a  decreased  to  aQ,  without  ever  producing 
large-scale  separation  of  the  boundary  layer. 

However,  the  interaction  with  the  inviscid  flow  pro¬ 
duced  the  distortions  in  the  fluctuating  pressure 
distribution  on  the  rear  half  of  the  airfoil  that 
are  evident  in  Fig.  It  should  be  mentioned  that 

these  quantitative  pressure  data  may  suffer  somewhat 
from  the  effects  of  tunnel  side-wall  boundary- layer 
contamination.  Nevertheless,  they  give  a  good  quali¬ 
tative  Indication  of  the  effects  of  mild  trailing- 
edge  separation  on  this  airfoil  and  on  the  VR-7 
airfoil . 


IV .  _  StaU  Onset :  Mild  Interact  ion 
Cene  r a  1  Ken  t  u  r es 

This  regime  has  special  practical  signifi¬ 
cance,  as  discussed  in  Sec.  II.  For  the  purposes  of 
mathem.it leal  modeling,  the  no-stall  and  stall-onset 
regimes  are  similar,  since  the  viscous  layers  remain 
relatively  thin.  However,  the  challenge  to  predic¬ 
tion  methods  is  greater,  because  the  magnitude  of 
the  viscous-inviscid  interaction  increases  rapidly 
with  small  increases  in  umax  in  the  stall-onset 
regime.  Furthermore,  depending  on  the  airfoil  geom¬ 
etry  and  M,„,  the  extent  of  flow  separation  may 
either  increase  or  decrease  with  increasing  reduced 
frequency,  whereas  viscous  effects  tend  to  be  sup¬ 
pressed  by  unsteady  effects  below  stall.  As  a  result 
of  these  factors,  the  stall-onset  regime  Is  a  narrow 
but  Important  one  that  bridges  the  no-stall  and 
light -stall  domains. 

I'.ase  5:  NASA  0012,  a  •=  9°  +  5°  sin  u»t,  k  3  0.20 

In  the  unsteady  experiments  at  M,„  c  0.30,  a 
thin  tongue  of  reversed  flow  akin  to  that  of  Case  A 
was  observed  on  the  rear  of  the  NACA  0012  airfoil 
for  W  <  -  13.9°.  For  k  <  0.05,  this  situa¬ 
tion  also  existed  for  *  1A°,  hut  when  the 

reduced  t requeue y  was  Increased  above  0. 10,  unusual 
separation- l ike  boundary-luyer  disturbances  origi¬ 
nated  in  the  upper  leading-edge  region  and  propagated 
downstream.  The  boundary- layer  instrumentation 
indicated  neither  reversed  flow  over  the  whole  air- 
toil  nor  complete  separation  in  the  sense  of  a  large- 
scale  breakdown  of  the  flow,  but  the  thickening  of 
the  viscous  layer  was  sufficient  to  cause  significant 
distortions  in  the  upper-surface  pressure  distribu- 
t  Ions  and  in  the  hysteresis  loops  of  C,  and  Cjq 
versus  a.  This  is  indicated  by  the  hold  curves 
In  Fig.  9 . 


of  the  figure  are  different  from  those  ot  the  pre¬ 
ceding  no-stall  case,  ap,M«  -  13”,  shown  in  fig.  M, 
but  the  effects  on  the  out-of-phase  and  second' 
harmonic  components  are  particularly  noteworthy. 

Case  6:  Vertol  VR-7,  a  *  10a  4  5^  sin  **t  , 

0.025  *  k  *  0.20 

In  contrast  to  the  NACA  0012  set  lion,  .1  well- 
defined  thin  layer  of  reversed  viscous  flow  devel¬ 
oped  in  the  trai ling-edge  region  in  the  si al 1 -onset 
and  light-stall  regimes  with  the  VR-7  airfoil.  K«*i 
Case  6,  the  chordwise  extent  of  this  reversed  flow 
was  strongly  dependent  on  reduced  frequency.  This 
is  illustrated  in  Fig.  11,  which  shows  the  loci  ot 
the  point  tw  r  0  on  the  airfoil  surface  as  .1  func¬ 
tion  of  time  for  three  values  of  k.  The  reversed- 
flow  region  is  clearly  delayed  and  suppressed  by 
increasing  unsteady  effects. 

The  upper-surface  pressure  distributions  are 
qualitatively  similar  over  the  experimental  range  of 
frequencies,  as  shown  in  Fig.  12.  However,  Fig.  13 
shows  that  the  extent  of  the  trai ling-edge  separa¬ 
tion  has  an  important  quantitative  effect,  and  only 
the  example  at  k  «  0.20  can  be  considered  truly 
representative  of  stall  onset.  For  k  <  0.10,  the 
results  exhibit  many  of  the  characteristics  ot  the 
next  level  in  the  hierarchy,  light  stall. 


V.  Light  Stall:  Strong  int e rac t ion 
General  Features 

This  regime  for  oscillating  airfoils  shares 
some  of  the  general  features  of  classical  static 
stall,  such  as  a  loss  of  lift  and  significant 
increases  in  drag  and  nose-down  pitching  moment  com¬ 
pared  with  the  theoretical  inviscid  values,  when  » 
exceeds  a  certain  value.  In  addition,  the  unsteady 
stall  behavior  is  characterized  by  large  phase  'a/s 
and  hystereses  in  the  separation  and  rent  t acluneni  ol 
the  viscous  f  l ow  and,  consequent  ly,  in  the  .111  loads. 
Also,  the  aerodynamic  damping  in  pitch,  given  bv 
♦Y^dn,  can  become  negative,  that  is,  unstable;  (hi-, 
tendency  is  strongest  in  the  light -stall  domain. 

Another  distinguishing  feature  oi  this  regime 
is  the  scale  of  the  interact  ion.  The  vertical 
extent  of  the  viscous  zone  tends  to  remain  on  tin- 
order  of  the  airfoil  thickness,  generally  Una.  than 
for  static  stall.  Consequent  ly ,  tills  class  ot 

airfoil  problems  should  be  within  t  )>• 
scope  of  zonal  methods  or  thin- layer  Navler-St  okes 
calculations  with  relatively  straightforward  tmbu 
lence  modeling. 

The  qualitative  behavior  of  light  stall  is 
known  to  he  esper i a  1 1 y  sens i 1 1 ve  t o  a i rl o l l  geom¬ 
etry,  reduced  frequency,  maximum  incidence,  and  Mach 
number;  also,  three-dimensional  effects  and  the  type 
of  motion  are  probably  important.  The  quant i l at i \ e 
behavior  Is  closely  related  to  the  boundary- lave* 
separation  character ist ics  (tor  example,  leading 
edge  versus  t rai l ing-edgo  separation!,  ami  to  tin- 
changes  in  this  separation  behavior  with  \nMN,  k, 
and  Mr,..  Therefore,  the  possible  test  cases  are 
numerous  and  varied.  For  the  sake  of  simplicity, 
however,  this  section  will  concentrate  on  the 
effects  of  reduced  frequency  on  two  airfoils:  one 
airfoil  with  trailing-edge  separation  character!;; 
tics  and  one  with  leading-edge  separation  behavior. 


Figure  10  shows  the  harmonic  components  of  the 
pressure  distribution  for  this  case.  All  tour  parts 


The  former  Is  represented  by  Case  6.  discussed 
earlier  and  presented  in  Figs.  11-13.  The  important 
point  here  is  the  suppression  of  the  amount  of 
trailing-edge  separation,  and  hence  the  magnitude  of 
the  viscous-inviscid  interaction,  with  increasing 
frequency.  As  a  result  of  this  unsteady  effect. 

Case  6  passes  from  light-stall  behavior  for  k  £  0.10 
to  stall-onset  conditions  for  k  >  0.15. 

Case  7:  NACA  0012,  a  -  10°  +  5°  sin  ut, 

0.025  US  0.20 

For  M„  <  0.2,  the  light-stall  behavior  of  the 
0012  airfoil  followed  the  trends  of  Case  6;  that  is, 
separation  was  suppressed  as  k  increased.  However, 
at  M„  ■  0.30  the  severity  of  the  stall  increased 
with  Increasing  k,  as  the  mild  trailing-edge  separa¬ 
tion  at  low  frequency  gave  way  to  the  separation- 
like  leading-edge  disturbances  that  were  noted  for 
Case  5  for  k  >  0.15. 

The  attendant  changes  in  the  lift  and  pitching 
moment  behavior  are  shown  in  Fig.  14,  which  provides 
a  striking  contrast  to  the  VR-7  data  in  Fig.  13. 

The  upper-surface  pressure  disturbances  are  given  in 
Fig.  15.  For  k  ■  0.025,  the  curves  generally 
resemble  those  of  Fig.  12,  but  the  results  at 
k  «  0.20  are  an  exaggeration  of  those  shown  in 
Fig.  9. 

The  maximum  lift  rose  monotonical ly  with  k 
for  both  the  0012  and  VR-7  airfoils,  as  shown  in 
Fig.  16.  However,  the  figure  shows  that  the  trends 
of  CM  are  quite  different  for  k  >  0.10,  as  the 

in  in 

stall  became  more  severe  for  the  0012  section.  In 
fact,  the  behavior  of  C^  versus  a  at  k  *  0.20 
and  the  near-collapse  of  the  leading-edge  suction, 
Fig.  15b,  are  characteristics  that  approach  those 
of  deep  dynamic  stall,  which  is  the  next  level  in 
the  hierarchy. 


VI.  Deep  Dynamic  Stall:  Viscous  Dominated 
Ccneral  Features 

The  time-dependent  stall  behavior  in  this  regime 
is  characterized  by  the  shedding  of  a  large  vortex¬ 
like  disturbance  from  the  leading-edge  region  and 
the  passage  of  this  vortex  over  the  upper  surface  of 
the  airfoil.1’  This  produces  values  of  Cj,,  CM, 

and  C|)  that  are  far  in  excess  of  their  static 
counterparts  when  a  is  increasing,  and  large  amounts 
of  hysteresis  occur  during  the  rest  of  the  cycle. 

The  scale  of  the  Interaction  zone  is  also  large,  the 
thickness  of  the  viscous  layer  is  of  the  order  of 
the  airfoil  chord  during  the  vortex-shedding  process. 
This  poses  potential  problems  in  grid  generation, 
spatial  resolution,  and  turbulence  modeling  for 
numerical  analyses. 

This  qualitative  picture  appears  to  hold  over 
a  wide  range  of  unsteady  flow  conditions,  airfoil 
geometries,  Reynolds  numbers,  and  Mach  numbers,  pro¬ 
vided  strong  shock  waves  do  not  develop.1’7  The 
quantitative  behavior  of  the  flow  depends  primarily 
on  the  time  history  of  the  angle  of  attack  for  the 
portion  of  the  cycle  when  a  exceeds  the  static 
stall  angle,  asa.  This  feature  and  other  details  of 
the  flow  behavior  will  become  evident  in  the  follow¬ 
ing  examples. 


Case  8:  NACA  0012,  q  -  10°  ±  10°  sin  wt ,  k  -  0.10 

The  right-hand  section  of  Fig.  2  shows  one  of 
the  examples  of  deep  dynamic  stall  for  the  NACA  0012 
airfoil.  The  contrast  with  the  static  behavior  and 
with  the  other  unsteady  stall  regimes  is  readily 
apparent. 

Figure  17  shows  the  unsteady  pressure  distribu¬ 
tion  on  the  upper  surface  for  this  case.  In  partic¬ 
ular,  the  vortex-shedding  phenomenon  manifests 
itself  in  the  curves  at  a  =  16.2“,  18.2°,  and 
19.4°;  and  the  leading-edge  suction  collapses 
abruptly  over  this  interval.  The  boundary  layer 
separated  initially  from  the  leading-edge  region  at 
this  frequency,  starting  at  a  -  16°,  or  wt  ;  40°. 
This  represented  a  change  from  abrupt  trailing-edge 
stall1  at  k  <  0.05. 

Case  9:  NACA  0012,  a  3  15°  +  5°  sin  uit, 

0.0Q5  <  k  <  0.20 

One  particular  reduced  frequency  of  this  set 
of  data,  k  ■  0.15,  produced  results  almost  identical 
to  those  of  Case  8  over  the  portion  of  the  cycle  for 
which  a  >  a  .  Figure  18  shows  the  time  histories 
of  a  ,  C^i  end  Cjj  and  the  corresponding  hysteresis 
loops  of  Cj^  and  C^j  versus  a  superimposed  on 
those  of  Case  8.  The  agreement  between  these  cases 
is  remarkable;  but  it  should  be  emphasized  that  the 
results  match  this  well  only  when  the  Lime  histories 
of  a  agree.1 

Figure  19  shows  the  growth  of  to  and  0,  , 

'max  ‘-min 

wii.li  reduced  frequency  for  the  NACA  0012,  along  witli 
the  Vertol  VK-7  results  of  the  following  case.  The 

extreme  values  of  Cm  ,  are  particularly  noto- 
nmin 

worthy,  as  they  are  considerably  larger  than  the 
light-stall  results  shown  In  Fig.  16. 

Case  10:  Vertol  VR-7,  a  =  15“  +  5"  sin  .t, 

0.025  iki  0.20 

At  lower  Mach  numbers,  the  VR-7  airfoil  exhib¬ 
ited  t rji 1 fng-edge  separation  behavior  throughout 
the  light-stall  and  deep-stall  regimes.  However, 
at  M.„  ~  0.30  an  interesting  change  occurred  in  tin 
boundary-layer  characteristics  from  trailing-edge 
separation  to  a  mixed  leading-  and  trailing-edge 
stall,  as  the  vortex  shedding  phenomenon  increased 
In  intensity  with  increasing  k.  The  lift  and 
moment  data  for  this  case  are  shown  in  Fig.  20  and 
the  boundary-layer  results  are  shown  in  Fig.  21. 

The  effects  of  vortex  shedding  first  appeared 
at  kg  0.05,  with  a  weak  vortex  originating  around 
midchord.  The  flow  reversal  that  preceded  this  vor¬ 
tex  formation  progressed  slowly  upstream  from  the 
trailing  edge,  as  indicated  in  Fig.  21.  However,  a 
conversion  to  leading-edge  stall  began  for  k  »  0.1') 
with  a  much  stronger  vortex  originating  at  x  ^  0.02. 
It  is  the  distortions  in  the  chordwise  pressure 
distribution  due  to  this  transient  vortex  that  pro¬ 
duce  the  large  values  of  the  pitching  moment  thal 
are  indicated  in  Figs.  19  and  20. 


VII.  Concluding  Remarks 

The  10  cases  described  in  the  previous  sections 
illustrate  most,  although  by  no  means  all,  of  ibo 
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panorama  of  viscous- inviscid  interactions  that  occur 
on  oscillating  airfoils  in  subsonic  flow.  The  spe¬ 
cific  flow  conditions  for  each  case  were  carefully 
selected  to  emphasize  the  predominant  characteristics 
of  each  of  the  four  flow  regimes,  while  minimizing 
the  effect  of  the  experimental  uncertainties  in  the 
appropriate  data.  In  other  words,  the  experimental 
results  were  presented  more  on  the  basis  of  document¬ 
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than  on  the  importance  of  the  absolute  values  of  the 
measurements.  Nevertheless,  the  data  are  thought  to 
be  of  sufficiently  high  quality  that  Including  the 
small  effects  of  the  wind-tunnel  walls  would  be 
warranted  in  any  analysis  of  the  test  cases. 

To  the  authors*  knowledge,  few  of  the  specific 
cases,  or  their  equivalents,  discussed  in  this  paper 
have  yet  been  calculated.  However,  the  unstalled 
cases  of  Sec.  Ill  would  appear  to  be  well  suited  to 
straightforward  extensions  of  present-day  zonal 
modeling  methods  for  steady  flows,  as  exemplified  by 
the  recent  collection  of  papers  in  Ref.  10.  Such 
methods  should  also  be  capable  of  treating  the 
stall-onset  cases  of  Sec.  IV,  although  the  details 
of  the  unsteady  flow  regions  with  flow  reversal  at 
the  bottom  of  the  boundary  layer  will  surely  require 
special  treatment.  The  thin- layer,  Reynolds-averaged 
Nav ler-Stokes  methods  (e.g..  Ref.  11)  are  also 
appropriate  for  this  regime,  as  well  as  for  the 
light-stall  regime.  The  extent  to  which  the  zonal 
methods  can  be  applied  to  the  light-stall  cases  is 
an  open  question.  However,  both  approaches  offer 
attractive  alternative  to  the  empirical  correlation 
methods  that  are  currently  used  for  engineering  pre¬ 
dictions,6*7  For  this  reason,  coupled  with  its 
practical  Importance,  the  light-stall  regime  is 
probably  the  one  that  most  warrants  concentrated 
research  efforts. 

The  greatest  challenge  to  theoretical  or  numer¬ 
ical  analysis  is  clearly  the  deep-stall  regime.  As 
mentioned  in  Sec.  VI,  the  viscous  region  is  very 
large,  and  it  is  a  dominant  feature  of  the  flow  field 
so  that  thin- layer  approxim.it ions  and  turbulence 
models  based  on  boundary- layer  length  scales  are 
probably  suspect.  New  approximate  or  semiempirical 
methods  that  take  advantage  of  the  predominance  of 
the  vortex-shedding  phenomena  might  be  feasible,  in 
addition  to  the  d iscrete-vortex  methods6*7  that  are 
sometimes  applied  to  bluff-body  flows,  provided  they 
are  not  limited  by  the  assumptions  that  have  to  be 
made  to  start  the  calculat ions.  The  long-range  need 
seems  to  be  far  Navi er-Stokes  calculations  with 
improved  turbulence  modeling. 
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Fig.  1  Airfoil  cross  sections. 
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Fig.  2  Four  regimes  of  unsteady  airloads  on  the  NACA  0012  airfoil;  solid  lines  denote  inct 
angle  of  attack,  dashed  lines  decreasing  a,  dotted  lines  static  data. 


STEADY:  U-V,  +  Vc  +  Vtcr 

UNSTEADY:  U  -  V,  +  Vc  ♦  V,  e>0  +  V,  a,  f  (x.O 

C-«  1  -  u2-^+ar(M4) 
p  at  » 


Fig.  3  Sketch  of  the  flow  due  to  thickness,  camber, 
and  angle  of  attack  on  an  oscillating  airfoil. 
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Fig.  5  Harmonic  components  of  unsteady  pressures  for  Case  1;  NACA  0012  airfoil,  a  -  3°  +  3° 


Fig.  12  Unsteady  pressures  for  Case  6;  VR-7  airfoil,  -i  *  10  +  V  sin 
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Fig.  13  Unsteady  airloads  for  Case  6 


Fig.  17  Unsteady  pressures  for  Case  8;  NACA  0012 
airfoil,  a  -  10°  +  10°  sin  u>t,  k  -  0.10. 


Fig.  18  Unsteady  airloads  on  the  NACA  0012  airfoil.  Cases  8  and  9 


Fig.  21  Loci  of  boundary- layer  flow  reversal  on 
the  VR-7  airfoil.  Case  10;  a  *  15°  +  5°  sin  »it. 


